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I. INTRODUCTION 



Important developments have been taking place in the growth of strained heterostructures based on Si and Ge^jSii-a; 
layers The particular case of Si layers pseudomorphically grown on relaxed Ge^^Sii-^; buffers is rather interesting 
because of the large mobilities reported, exceeding the corresponding bulk Si values |^,||. It has been shown that in 
this type of double heterostructures (DHS) the otherwise degenerate six Ag valleys of bu k Si are shifted in energy. If 
the DHS is grown along a high symmetry direction (say, the < 001 > direction), the two valleys along this direction 
are shifted downwards in energy, while the other four valleys are shifted upwards. As a result, we are led to two 
A2 valleys with a bottom energy below the corresponding one for the bulk valleys and four A4 valleys with a higher 
bottom energy. The energy shift between the A2 and A4 valleys is empirically estimated as AE — 0.6 x eV. In this 
kind of DHS, a quantum- well like band alignment in the Si channel was unexpectedly found in 1985 which is 

due to tensile strains in Si and ensures high-mobility n channel for the doped structure. High mobilities are explained 
as a consequence of the low effective mass of the carriers in A2 valleys and also because the scattering efficiency of 
the carriers (by phonons and impurities) is reduced. The above mentioned facts are of great importance for device 
performance and this kind of structures seems to have a prominent future for applications in micro and optoelectronics 

i- 

Monte Carlo calculations of high electric field mobilities and drift velocities for electronic transport along the Si 
channel of Si/Ge^^Sii-^^ were reported by various authors [^-[sj and also applied in the study of modulation-doped 
field-effect transistor (MODFET) structures ||]. More recently, hole transport parameters have been analyzed for this 
type of systems 0. While the electronic transport properties of low field, low temperature Si/GeSi heterostructures 
are more or less well understood the situation is not the same for high electric fields. For high field transport, 
the effects of size quantization are usually negligible and we can work within three dimensional models for the valley 
structure. In Ref. [Q, size quantization was introduced in a 10 nm Si channel, but the authors had to include up to 
six subbands in the A2 valleys. Strain effects are introduced through the splitting energy AE between the A2 and 
A4 valleys, which ranges from 0.1 eV to 0.4 eV. Such relatively large values of the energy shift effectively reduce the 
intervalley phonon scatterings between the valleys (in comparison with the unstrained Si). 

In the present paper we report results for high electric field electronic transport along the Si channel of a Si/GeSi 
DHS. We model the system in close analogy with the above mentioned works, taking into account the energy shift 
between the A2 and A4 valleys and neglecting size quantization. Rather than Monte Carlo simulations, we apply 
an analytical Fokker-Planck approach (FPA), which treats transport as an energy-diffusion process in energy space. 
This approach was proposed a long time ago [ ^Jl3| as an alternative for the consideration of high-field transport in 
semiconductors. The theory bears a semiclassical nature and has been recently revisited, the general formalism being 
discussed with further details ||l^-[l6|]. As a test of the FPA in the case of a well known semiconductor, the present 
authors applied the theory to bulk Si, where both experimental data and Monte Carlo simulations for the high field 
drift velocity were successfully reproduced |]l7| . The current calculations can also be considered as a way of testing 
the FPA in a somewhat different system, which is of present day interest in high technology applications. 

The FPA is applicable when the energy exchanges between the carriers and the surrounding medium (crystal lattice 
-|- external field) can be assumed quasicontinuous. This latter condition apparently invalidates the method for highly 
inelastic scattering mechanisms (as is the case of carrier scattering by optic and intervalley phonons). However, if the 
carrier energy is large enough, the exchanged energy becomes certainly low compared with the energy of the carriers, 
and this is the case of high field transport. Hence, the FPA is assumed to be valid if the condition Eav » fiio 
is fulfilled, where Eav is the average carrier energy and hio the exchanged energy (say, the phonon energy). The 
FPA has the advantage of being analytical, and, whenever it can be applied, saves computational time and allows a 
more transparent physical interpretation. Of course, it cannot compete in accuracy with Monte Carlo simulations. 
However, we have found that the FPA leads to comparatively good results even in cases where several scattering 
mechanisms should be taken into account ||l^. In contrast with Monte Carlo simulations, the analytical approach 
involves less realistic models for the analysis of a concrete semiconductor. More details about the applied model and 
the fundamental theory are given in the next sections. 

The paper is organized as follows. In Sec. II we briefiy summarize general theoretical aspects of our work, in Sec. 
HI details of the calculations for the Si/GeSi DHS are presented, in Sec. IV the results of our calculations are shown 
and comparison with previous works is made. 



II. GENERAL THEORY 



The FPA considers transport in the spirit of a drift-diffusion process in energy space. A certain distribution function 
(DF), f{E, t), is defined which depends on the carrier energy E and the time t, such that f{E, t)N{E) gives the number 
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of carriers at time t having their energies in the interval [E, E + dE]^ while the function N{E) represents the density 
of states (DOS). The DF obeys the equation ||-|6| : 



where 



JiE, t) - W(E)N{E)f{E, i) - ^ [D{E)N{E)f{E, t)] . (2) 

In Eq. (H) W{E) represents a certain "drift velocity" of the carriers in energy space and in fact gives the rate of 
energy exchange of the carriers with the surrounding medium, while D{E) is a kind of diffusion coefficient. Equation 
(|l|) has the form of a continuity equation for the carrier "motion" in energy space and leads to the conservation of 
the carrier "flux" . For a many valley system it should be applied to each of the valleys separately and it is worth to 
note that Eq. (|^) does not describe intervalley coupHngs. The quantity J{E,t), given by Eq. (||), is thus the carrier 
current density in energy space. Under steady-state conditions the DF is time independent. Moreover, we must have 
J{E) = 0. Hence, we are led to 

A [DiE)N{E)f{E)] - W{E)NiE)fiE). (3) 

We assume the existence of a phonon bath in equilibrium at the temperature T. The carriers interact with the 
phonons and the applied dc electric field F. We suppose, as an approximation, that the continuous exchange of phonons 
between the carriers and the phonon bath does not affect the thermal equilibrium of the latter. The coefficients W{E) 
and D{E) are split as follows 

D{E)^DF{E)+Dph{E) , W{E) = WFiE) + Wph{E). (4) 

The label "F" {''ph") denotes the electric field (phonon) contribution to these coefficients. We are obviously neglecting 
the intracoUisional field effect. The explicit form of these coefficients is obtained below. 
Equation (^) has the simple solution 



(5) 



where Dph{E) was neglected. This approximation is very well fulfilled in all the cases of interest for us p|,|3. The 
practical usefulness of the FPA lies in the possibility of the analytical performance of the integral in Eq. (m). 

Let us apply the above formalism to electrons in the A2 energy valleys of strained Si of a Si/Gea;Sii_a; DHS. The A2 
valleys are ellipsoids of revolution in fc-space with their revolution axis along the < 001 > crystallographic direction, 
coincident with the growth direction of the DHS. To be specific, let us take F = (0, F, 0), where the x, y, z coordinates 
are taken along < 100 >, < 010 > and < 001 > respectively. The energy dispersion relation is given in the form 
e = e{p). In order to take into account non-parabolicity of the band structure, we assume 

7(6) = eil + ae)^:^ + ^^^. (6) 
Znit 2mi Zmo 

In Eq. (|^) rut (mi) is the transverse (longitudinal) effective mass of the bulk Si conduction band electrons. We suppose 
the masses are not changed by the strains in the Si layer: mt/m^ — 0.19 and mi/mQ — 0.916. Nonparabolicity is 
estimated also the same as in bulk Si: a — 0.5 (eV)^^. The Herring- Vogt transformation was applied in writing the 
right-hand side (RHS) of Eq. (^, where pt = {mt/moy/'^p^, pi = {mi/moy/'^p'^ and p*"^ = p^"^ +P;*^- 

The explicit expressions for Dp{E) and Wph{E) can be directly taken from Refs. |p"^-jl5[|. Hence, for Dp{E) we 
have 

DAE)^'^4^T{eUe)/{i{e))\ (7) 
3mt 



where 7'(e) denotes the first derivative of the function (defined in Eq. (|6|)). For r(e) we shall consider one "5" 
intervalley phonon responsible for transitions between the equivalent A2 valleys. Then the relaxation time reads as 
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+ K(T) + - hLOg)\l + 2a(e - nwg)\e{e - hug) , (8) 

where 9{e) is the step-function, 

ng{T) = [eMf>^g/kBT)-ir\ (9) 

and 

3/2 n2 

Cg = -J (10) 

\/2Tiph?UJg 

In Eq. (^o|) p is the semiconductor mass density, ujg and are the phonon frequency and deformation-potential 
(DP) constant respectively for intervalley phonons of type "g". As a simplified model for our calculations, we have 
considered just one "(7" phonon with an energy fiiOg = 0.031 eV and DP coupling constant Dg = 11 x 10® eV/cm, 
which approximately corresponds to the average of the three "g" phonons of bulk Si reported in Table VI of Ref . . 
Other bulk Si parameters are also taken from the same reference. 
For ellipsoidal A2-valleys of the Si layer CB, the DOS is given by 



N{E) - -j^J^^Eil + aE)\l + 2aE\, (11) 

where nonparabolicity is taken into account and nid = \fm^nvi. In calculating Wph,{E\ we use the expression 

Wph[E) = hu, [l/rafc.(p) - l/re™(p)] , (12) 

where "a&s" ("em") denotes phonon absorption (emission) by the electron. We must remember that 

l/r(e) - l/T„bs(e) + l/T-e„,(e), (13) 

the "afts" ("em") term in Eq. ([l3| ) corresponds to the first (second) term at the RHS of Eq. (||) for the case of 
intervalley "g" phonons. Hence, using Eq. (^) and latter expressions for Tabs and Tem, we can find the explicit form 

Intravalley optical phonons do not contribute to transition rates because the corresponding transitions are forbidden 
by the selection rules. The contribution of intravalley acoustic phonons will be ignored in the present work. For high 
temperatures and high carrier energies they should provide a weak contribution to transport parameters and we shall 
limit ourselves to this case. 

As was remarked above, intervalley phonon scattering, strictly speaking, is beyond the scope of the Fokker-Planck 
equation in its standard form [Eq. (^]. After the absorption or the emission of an intervalley phonon, the electron 
is dropped away from one valley into another one, and such abrupt transitions are not explicitly assumed in Eq. (|^), 
where the electron number is conserved within each valley. However, when we have equivalent valleys (as is the case 
of A2 valleys), we can still apply the Fokker-Planck equation in its standard form; electrons dropped away from one 
valley fall into an equivalent one and vice versa, all happens as the electrons always remained within the given valley. 
Of course, if the valleys are not equivalent this argument no longer applies. This is the case of electron transitions 
induced by the "/" intervalley phonons between valleys A2 and A4, which are not equivalent in the strained Si layer 
of the Si/GeSi DHS. For such intervalley electronic transitions we must consider a different approach. 

In order to estimate the effect of electron intervalley transitions due to the interaction with "/" intervalley phonons 
we must realize that this effect is actually weak and should introduce just small changes into our final results. As for 
the case of "g" intervalley phonons, we assume just one type of "/" phonons with energy fiLOf = 0.042 eV and DP 
coupling constant Df = 4 x 10® eV/cm (we again estimated these parameters as averages from the three corresponding 
"/" phonons reported in Table VI of Ref. [|8[). In all cases we considered zero-order intervalley phonons in our model 
(a more realistic model must consider first order intervalley phonons in the way discussed in Ref. [p9|). Electron 
transitions due to the interaction with "/" intervalley phonons from A2 valleys up to A4 valleys is formally switched 
on at i = and the DF in each of the A2 valleys evolves in time by the approximate law ^ = —J/tj , where tj 
and the corresponding expressions for nf{T) and C/ are obtained from Eqs.(g), (^ and ([l0| ) by means of the formal 
substitution: g f and e ^ e — AE. 

After a time r has elapsed, we obtain the result 

f{E,T)=f{E)exp[~T/TfiE)], (14) 
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where f{E) represents the DF without taking the A2-A4 intervalley scattering into account. Of course, this is just a 
rough estimation of the effect and we are actually neglecting the inverse intervalley process (from A4 to A2 valleys). 
Due to the weakness of the effect, we shall assume that this procedure is satisfactory. For the parameter r in Eq. 
(|r^), we make the reasonable estimate 

1 



rrithhjf 



2n.(T) + l 



(15) 



III. CALCULATION OF THE DISTRIBUTION FUNCTION 

From Eq. (^, with the explicit application of Eqs. (0) and (12), we are led to the DF for A2 valleys without the 
consideration of the A2-A4 intervalley scattering process. The final result is 

J,{E)=E^^{l + aE)^^eMPtEP,{E,T)), j = 1, 2, (16) 

where Aj and Bj are parameters (dependent on T and F) and Pj{E,T) is a polynomial in E. This structure is far 
from the Maxwellian one. The DF describes a stationary non-equilibrium configuration where an electron temperature 
Te cannot be defined. In Eq. (|l^) we shall measure all energies in units of huig and 

inif'm'iDt 

For £' < 1 we obtain 

= Pt{ng{T)f{l + Eo){l + 2E^f 

B^^Pt{ng{T)f{E^-l){Eo-l/2f 

4 

P^{E,T)^{ng{T)fY.''-^^'^ (18) 

1=0 

where the an are coefficients explicitly dependent on Eq — fkOga. For E > 1 the following result is obtained 

A2 - Pt [{ng{T))^2{l + 8E^) - {2ng{T) + 1){AEI - SEi + 5Eo - 1)] 
B2 = Pt [{ng{T))^2{8E^ + 1) + i2ng{T) + 1)(4S3 ^ 3^2 ^ ^ 

4 

P2(i^,T) = ^a,2£^\ (19) 

1=0 

where the coefficients 0^2 are now dependent on T through ng{T). All the coefficients Uij are given in Table I. 

In the interval E' > 1 we must distinguish three subintervals: 1 < < AE — Ef, AE — Ef < E < AE + Ef and 
E > AE + Ef, where Ef = huf/hujg. In the latter two subintervals, the DF shall be given by 

f,{E,T)^f2{E)exp[~T/Tf{E)] , J =3,4, (20) 

in correspondence with our approach to intervalley transitions between A2 and A4 valleys. 

Once we have determined the DF, we can calculate the electron average energy by the expression 

Eav = J {EJ{E)N{E)}dE/ J {f{E)N{E)}dE. (21) 

In Eq. (|2^) we just consider A2 valleys, i.e., we assume the A4 valleys essentially deprived of carriers, an approximation 
that seems acceptable if the applied electric field is not extremely large. Moreover, in Eq. (|2^) we should be careful 
in partitioning the integration into the four energy subintervals mentioned above. 
Another important quantity is the drift velocity Vd given by 



2eF f \i{E)Tg{E) 

Vd ^ 



3mtJ { {l'{E)y 



dE 



N{E)jdE/ I {f{E)N{E)}dE. (22) 



In Eq. (|2^) we must partition the integration in the same way as in Eq. (pl|). We just considered the drift velocity 
from electrons in the A2 valleys. 
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IV. DISCUSSION OF RESULTS 



By direct application of the theory developed in the foregoing sections, we have made numerical calculations of 
the average electron energy Eav and the electron drift velocity Vd as functions of the temperature T and dc electric 
field F. Our results are essentially valid for high temperatures as far as intravalley acoustic phonon scattering was 
ignored. As it was remarked in Sec. I, the FPA is applicable for high electric fields, when the condition Eav >> fno 
fulfils In the examined Si/GeSi DHS, we take into account the relatively large energy shift between valleys 

A2 and A4 and consider that just the A2 valley is substantially populated by carriers. Hence, intervalley transitions 
between valleys A2-A4 are assumed to be weak processes and treated within a relatively coarse approximation. As it 
was discussed in the previous sections, such transitions are actually beyond the scope of standard FPA. All numerical 
parameters used in computations were taken from Ref. fl^ and also shown in Sec. III. 

In Fig. 1 we show our results for the average electron energy Eav as a function of electric field F for T — 300 K. 
Two values of AE were considered: 0.4 and 0.1 eV. As expected, the results show a rather weak dependence on AE. 
The involved carrier energies are large enough, thus ensuring the applicability of the FPA. After comparison of Fig. 1 
with Fig. 1(a) of Ref. Q, we can see that a reasonable agreement was actually achieved. It is important to notice 
that just one fitting parameter was applied in the f3t of Eq. (p^), which could be related to the overlapping integral 
describing the electron-phonon scattering probabilities. This overlapping integral was assumed unity in the general 
formulas of Sec. II, but actually it differs from unity when intervalley phonons are present. The linear behavior oi Eav 
for < 5 kV/cm should not be realistic because the FPA is not valid in the low-field regime as discussed exhaustively 
in previous works. ||l^[l^ 

In Fig. 2 we present the drift velocity Vd as a function of the electric field for T = 300 K. The continuous curve 
represents our calculations, while the dots were taken from Fig. 2(a) of Ref. jl^ . We just considered the case AE = 0.2 
eV. As a matter of fact, the curves for different values of AE are almost coincident (in close agreement with the results 
shown in Ref. |15|| ) as can be expected from general physical grounds. For lower electric fields our results deviate from 
those of Ref. [ [15[ , a reasonable result taking into account that the FPA is valid for high electric fields. However, for 
very high electric fields {F > 80 kV/cm) we again notice an increasing deviation from the Monte Carlo results of Ref. 
[ p^ . The FPA is not able to describe the saturation value of Vd in the very high electric field region. On the contrary, 
it is obtained a decreasing behavior of the drift velocity at a rate that becomes much higher as higher is the field. In 
our present treatment (as well as in that of Ref. |jl^) negative differential mobilities are out of question as far as the 
contributions of electrons from the A4 valleys are not considered. As a final remark, we should notice a very good 
agreement with the results of Ref. jlH) for a wide interval of electric fields. 

In Fig. 3 we show the same plots as those of Fig. 2, but now for T = 77 K. From a simple glance at Fig. 3, it is 
obvious that for such a low temperature the agreement between the FPA and the Monte Carlo results of Ref. jl^] 
(Fig. 2(b)) is much worse. As it was discussed before, we have ignored intravalley acoustic phonons. Even though 
this is not a limitation of the FPA itself, it is a difficult task to include them in the calculations. Furthermore, at low 
temperatures, quantum effects should become relevant and a quantum approach should be required. Hence, we stress 
that the FPA results are reliable just for high temperatures. 

In this work we have compared our results with those from Monte Carlo simulations. ||l^ Other possible comparisons 
could be done. For instance, we could compare with Fig. 3 of Ref. |^ or Ref. 1^. But all these results are in accordance 
with themselves and nothing essential would be added. Our results are also comparable to those of Ref. where 
size quantization was examined. As it was said before, size quantization is of relevance when the carrier energies 
are low enough (this is the case of low field transport as discussed in Ref. [0), but for the high carrier energies 
involved in high field transport, size quantization becomes irrelevant. In the revised literature we have not found 
available experimental data for this kind of system. As a general remark, we conclude that the FPA leads to results 
in acceptable agreement with those of Monte Carlo simulations in those intervals of temperature and electric fields 
where it is supposed to be valid. This conclusion is true in spite of the more or less coarse simplifications we have 
to face within the limits of this method. However, if wider intervals of electric fields or more accurate treatment is 
required, Monte Carlo simulations or some other equivalent numerical procedure should be necessary. 
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FIGURES 



FIG. 1. Average electron energy ( in eV) as a function of the electric field for T = 300 K. Two values of AE are 
examined: 0.4 and 0.1 eV. 

FIG. 2. Drift velocity as a function of the electric field for T = 300 K. Our results are represented by the continuous 
while the dots were taken from Monte Carlo simulations of Ref. We have set Ai? = 0.2 eV. 

FIG. 3. Same as in Fig. 2 for T = 77 K. The dots are Monte Carlo results from Ref. As one can see the 

agreement is now much worse. 



TABLE I. Coefficients aif. i = 0, ■ ■ ■ 4 , j = 1, 2 



4£o(4£;g + 8Sg + ll£o + 3) + 1 


8Eo{8E'^ + 3) + {2ng + l)(-16£;,t + 32£;g - 44Sg + 12Eo - 1) 


AEo{9,E^ + 12El + TEo + 1) 


8E'int,{8E'o +7) + {2ng + l)(8£o* - 12£o + 7Eq - 1) 


8£g(2So + l)^ 


64S,tn^ - 8Ef,{2ng + l){-iEf, - AEa + 1) 


8£g(2£;o + 1) 


32S(tn^ + 8Et{2ng + l)(2Eo - 1) 


16^0 /5 


-l&E^{2ng + l)/5 
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